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Abstract: 

We discuss the effect of 77-77' mixing and of axial current anomalies on various decay pro- 
cesses. We describe the 27 decays of 77 and 77', the radiative decays of low-lying vector 
mesons as well as the decays 77 — ► 3n, J/ip — > 777 and J/ip — > 77*7 in one unified frame- 
work, in which they are related to the electromagnetic and strong axial current anomalies, 
respectively. We also discuss briefly the enhancement of rj' in D s decays and suggest a pos- 
sible relation to gluon-mediated processes. The importance of the decays of the glueball 
candidates / (1500) and / (1590) into r/r/' channels is also stressed. 
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1. In this paper we study the effect of 77-77' mixing and of the axial current anomaly on 
various processes. We show that a consistent picture arises for the radiative decays rj —>■ 27 
and 77' — ► 27, for radiative decays of type V — > 777, rj' — > of the lowest-lying vector 
mesons V, as well as for 77 — > 37r and J/ip — > r)(r]')j. Whereas for all the other decays we 
follow closely the approach of Ref . || , our description of the vector meson decays is based 
on their relation to the QED triangle anomaly In Sec. 2 we thus first test our approach for 
the simpler cases of the decays V — > ^(K)^. In Sec. 3 we then establish notations for the 
77-77' system and determine the decay constants /o and f$ as functions of the mixing angle 
from the 27 decays of 77 and 77'. Section 4 is devoted to the comparison of our expectations 
for the widths of V — > 777, 77' — > V7 with experiment, whereas Sec. 5 deals with the decay 
77 — ► 37r, which probes the light quark content of 77. In Sec. 6 we then investigate the 
gluon content of 77 and 77', which determines the decays J ftp —>■ 777 and J ftp — > 7/7. It 
is in this channel that our approach yields a dependence of the involved matrix elements 
on the mixing angle that is significantly different from other models. We briefly comment 
on this in relation to the approach of Veneziano et al., Ref. 0, and close in Sec. 7 with 
some remarks on the important role of anomalies in other gluon-enriched channels like 
D s — > T](j]')X and in glueball searches. 

2. The radiative decays of the lowest-lying vector- meson nonet are traditionally described 
in terms of magnetic moments of quarks, see e.g. 0, or of unknown couplings related 
by SUf(3) symmetry [f|. In Refs. [§, § the width T(p — ► nj) was related via SUf(3) 
arguments and vector- meson dominance to the radiative width of tt. In this paper we 
relate the radiative decays of light vector mesons to light pseudoscalars, V — > P7 and 
P — > V7, directly to the anomaly of the AVV triangle diagram, where A stands for an 
axial-vector and V for a vector current. Our approach[| both includes SUp(3) breaking 
effects and fixes the vertex couplings gyp 1 as defined belowQ. 

Let us start by considering the correlation function 

ijdhe*** (n( qi + q 2 )\ T J™(x)^=°«(0) |0) = i^ =0(1) ( gi 2 , qj), (1) 

with the currents 

2 11 1 

J ™ = 3 ^ lflU ~ 3 ^ lfld ~ 3 ^ 7/ " S ' Jl ^ =01 = 71 ^ lllU ± ^ 7m ^ ' ^ 

The values of F n j(0, 0) are fixed by the QED triangle anomaly as 

^-(0,0) = ^, ^(0,0) = ^, (3) 

where the pion leptonic decay constant f n is defined as 

( I u>y^ 5 d I 7r~ > = if v pp, (4) 
X A similar approach was proposed in Ref. [Q. 

2 The behaviour of the anomalous contribution has been studied quite extensively, in particular for 
large space-like momenta, and found to fall off asymptotically like 1/ — q 2 || . The results of this section 
show however that in the physical sector, and for momenta up to the square of the <p mass, the full value 
of the anomaly gives sensible results. It was important for us to check this fact in the p, u>, n, K sector, 
independently of further difficulties associated with the strong anomaly. 
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so that f w = 0.132 GeV from n — ► /xty For the analogously defined kaon decay constant 
one finds fx = 0.160 GeV. Using their analytic properties, we can express these form 
factors by a dispersion relation in the momentum of the isospin current^: 

1 °° 

F WiJ=1(0 )(0,0) = - / ds-lmF nJ=m {s,0). (5) 

71 J S 

4(9)m2 

Saturating by the lowest-lying resonances we obtain 

F nJ=1 (0, 0) = A ^ 7 + . . . , F ff ,j =0 (0, 0) = ^ + ... (6) 

Here the dots stand for higher resonances and multi-particle contributions to the correlation 
function. In the following we assume vector meson dominance and thus neglect these 
contributions. Although this assumption may be criticized (and actually in Ref. || the 
continuum contributions were modelled by the perturbative spectral function above some 
threshold, but without taking into account the motion of the light quarks in the pion), it 
yields a description of the data that is good enough for our purposes^. 

The fy in (|5|) are the vector mesons' leptonic decay constants defined by 

(0| J^\V(p,X)) = m v f v e^(p). (7) 

A denotes the helicity state of the meson. For ideal mixing the relevant currents are 
J? = J* =1 , J 1 ^ = J^ =0 and = —s^^s. In the following we include the deviation from 
ideal mixing by taking into account a mixing-angle 9y = 40.3° for <fi-u mixing^. The fy 
can be determined from the experimental decay rates || via 

T(V -> e+O = cyira 2 (8) 

my 

with cy = {2/3, 2/9 sin 2 Q v , 2/9 cos 2 9y} for V = {p°, u, <fi}. The experimental results are 
f p0 = (216±5)MeV, f u = (174 ± 3) MeV, / = (254 ± 3) MeV. (9) 
The charged mesons decay constants can be obtained from r — > Vv T via 

F(r- - V«,) = !Wr U - 4V (l + 2 4) , (10) 

lo7r V J \ J 

where Vij is the appropriate CKM matrix element. We find 

f p ± = (195 ± 7) MeV, f K «± = (226 ± 28) MeV. (11) 

In view of the discrepancies between the measurements of the p decay constant from the 
charged and neutral sectors, we will use the average value f p = (205 ±17) MeV in the 
following. 

3 Note that the dispersion relation needs no subtraction, and that the value ^^(0,0) is unambiguously 
fixed by the anomaly. 

4 We are interested here in understanding how the gluon anomaly affects decays involving r\ and rf and 
not so much in a detailed fit of radiative vector meson decays. 

5 This value follows from the standard SUf(3) breaking analysis; the ideal mixing angle is given by 
tan fly = l/y2, corresponding to 6y = 35.3°. 
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P V 



7T ± p 1 * 1 



7T° p° 



7r o; 

K° K°* 
K ± x±* 



g VPl (th.) 



3m u sin 6y + \/2 cos Oy 



(0.73 ±0.07) GeV" 1 
(0.73 ±0.07) GeV" 1 

*UK* -P-58±0.04)GeV- 
^!" 9 (V2sin0 y - cosfly) = (0.20 ± 0.01) GeV" 1 
mK * - (1.25 ±0.15) GeV" 1 

(0.62 ±0.08) GeV" 1 



5fyp 7 (exp.) 
(0.74 ±0.04) GeV" 1 
(0.98 ±0.13) GeV" 1 
(2.33 ±0.07) GeV" 1 
(0.14 ±0.01) GeV" 1 
(1.25 ±0.05) GeV" 1 
(0.84 ±0.04) GeV" 1 



Table 1: Theoretical and experimental values of the on-shell V-P electromagnetic vertex 
couplings defined in Eq. (|T2"D. For gyp-y (th.) we give only the experimental errors coming 
from the decay constants fpy. 

Finally, we introduce vertex couplings gvp^y, which are just the on-shell V-P electro- 
magnetic form factors: 



PM I Jf" I V( Pv , A) ) 



{pv-Pp) 2 =0 



v p (*)(A)<7 

'dVP-y £[iupaPpPv € V 



'121 



The amplitude of the decay P — > V7 or V — > Pj, depending on kinematics, is then 
obtained by contracting with the polarization vector of the photon and multiplying by 
yivra. The decay rates read 



a 



T(P^V 1 ) = -g VP . 



m. 



m 



v 



mp 



Y{V -> Pi) 



a 
24 



9vp-) 



m 



v 



nil 



(13) 



In Table [I] we give both the theoretical and experimental values of the couplings gvp-y 
for decays involving pions or kaons. The agreement between theoretical and experimental 
values shows that the ground states indeed dominate the spectral functions in (|5|), with 
corrections of order 10%. There are, however, three channels where the deviation from the 
experimental value is larger. The couplings in the p-n channels should be equal, those 
in the K*-K channel differ by a factor of two, which is not quite supported by the data. 
Since the couplings are related by Clebsch-Gordan coefficients, we see no possibility to 
reconcile the data with our predictions and leave this "anomaly" as open question to the 
experimentalists. For <fi — >• 717 the deviation is more than 10%. This decay, however, is 
in our approach completely due to (fi-u non-ideal mixing and thus strongly suppressed 
compared to the other channels. It may also be influenced by other, usually negligible 
mechanisms, in particular p-oj mixing, which we have neglected here, so that we consider 
the agreement with the experimental coupling as still satisfactory. 
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3. Having gained some control over the radiative decays of vectors into tt and K, we now 
turn to the rj-rj' system, which is our central point of interest in this study Numerous 
possibilities, of which we can only quote a few JTOj , have been suggested to describe this 
system, with or without explicit mixing with extra "glueball" states. The discussion of all 
these approaches goes far beyond the scope of this letter, and we will content ourselves with 
shortly commenting in Sec. 6 on the elegant approach by Veneziano et al., Ref. [[|. For 
the time being, however, we focus on the more simple and pragmatic approach proposed 
in Ref. Following the Particle Data Group conventions, which differ from the ones 
originally used in Ref. 0, we reproduce here the key formulas in the current notation. In 
terms of the mixing angle 9 and their singlet and octet components rjo and r]s, respectively, 
the physical states are decomposed as 

\v) = \Vs) cos 9 — 1 770) sin^, |r/) = \r] 8 ) sin 9 + |r/ ) cos6>. (14) 

Traditional SUf(3) based analyses suggest either 9 ~ —10° or 9 ~ —23° from the quadratic 
and linear version, respectively, of the Gell-Mann-Okubo mass formula, whereas we leave 9 
as an open parameter, at least for the time being. Analogously to the pion decay constant 
f n , Eq. (f|), we define the decay constants^] / 8 and f as coupling of r/ 8 and r/ , respectively, 
to the divergence of the relevant axial- vector currents (where f$ = f n 7^ fo in the SUf(3) 
limit): 



d^Ag = -j= (m u uij 5 u + m d dij 5 d - 2m s si^ [ 

2 13 cy. 

d^Aft = -j= \m u ui^u + m d di^ b d + m s si^s} + -j= -~ G^ V G A ^ V ', (15) 



where G A is the gluonic field-strength tensor and G A = 7;€ f j, l/pr jG Apcr its dual. Defining 
as in Ref. [Ij the interpolating fields of 77 and v( as linear combinations of the axial-vector 
divergences, we thus have 

(0|<9Xh> =mlf 8 cos9, (0|dXh> = -m 2 v f sm9, 
(0\d lt A£\r/)=m*,f a sme, (0 | \ V f ) = m 2 v ,f co S 9. (16) 

Our essential assumption is that the u and d quark masses can be neglected in (|15"D. 
This yields the following simple expressions for the matrix element of the strong anomaly 
over the vacuum and //(?/) : 

(0\i?±GG\ri) = y|mj (/ 8 cos 9 - V2f sin 0) , 
(0\^GG\r)') = (f 8 sm9 + V2f cos9 



(17) 



6 Note that by virtue of the "hard" non-conservation of Aq, fa is scale-dependent as described in detail 
in the first reference of [0; we will assume here that the momentum dependence can be neglected. 
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Figure 1: fo/fn and fs/fn, determined from rj — > 27 and rj' — > 27, as functions of the 
mixing angle. 



We thus have three parameters to fix: fo, fs and 9. Two of them can be determined from 
the radiative widths of rj and rj', which are given by: 



3 

mz 



9 /cos6> 2v / 2sin6>\ 
96vr 3 V f 8 fo J 



m^, 2 (sm9 , 2 a/2 cos 6> 



2 



r C'-^ = 9^ a2 (lT + ^J' (19) 

As it was done in Ref. |T|], we solve these relations for the decay constants and thus get fo 
and fs as functions of the mixing angle. As experimental input we use || 

r(?7 -> 2 7 ) = (0.51 ±0.026) keV, F(rf -> 2 7 ) = (4.53 ± 0.59) keV, (20) 

where according to the suggestion in the full listings part of the Review of Particle Prop- 
erties (Ref. |§, p. 1451) we have only retained the more recent data. In Fig. [L], fo/fn and 
/s/Ztt are plotted as functions of the mixing angle (solid lines) together with their exper- 
imental error (dashed lines). We find thus fs/ fn > 1 for 9 less than pa —15° as expected 
from chiral perturbation theory. The precise value |TTJ fs/ fn — 1-25, however, corresponds 
to a rather largeQ mixing 9 = —21.3°. 

Having thus fixed fo and /§, we continue in the next section with the electromagnetic 
properties of 77 and rj', i.e. the vector meson decays, consider then in Sec. 5 their light quark 
matrix elements, to finish in Sec. 6 with their glue content. 

4. We are now in a position to deal with the radiative vector meson decays involving rj 
or v( . In Table |2| we list the decay channels together with the theoretical formulas for the 
couplings gyp-y and their experimental values. In Fig. we then plot the theoretical and 
experimental values as functions of the pseudoscalar mixing angle. Since in our approach 
continuum contributions are neglected, the couplings should be slightly overestimated by 

7 In fact, Ref. [H supplements the more traditional approaches to the 27 decays of rj and 77' by taking 
into account continuum contributions, which just cancel the effect of the large SUp(3) breaking ratio 
fs/ fit = 1-25. The result is then 9 = — (17±2)° with the experimental input data of Eq. (p0|). As we shall 
see in the following, this value is in excellent agreement with the results we find from the investigation of 
other decay channels. 
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Figure 2: Theoretical and experimental values of the V-rj, rf electromagnetic couplings as 
functions of the pseudoscalar mixing angle. Solid lines: theoretical predictions according 
to Table 0. Short-dashed lines: errors of the theoretical couplings coming from the exper- 
imental errors of fo and /§. Dashed-dotted lines: experimental couplings including errors 
(for g^'-y, plotted in (f), there exists only an experimental upper bound). 
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T]' <f) 



., ray 
g VPl m units - — - 
/y7T 



4/ 8 



tv cos y — 



|3 1 , V3 . , 

«-r cos ^ + 7T sm ^ 



3 1 
v/3 



— sin# 
/o 



4/ 8 



2 cos 6V — sin 9y) ~f=~ sm @v 

2v2/o 



/. 

cos# 

cos8 . n sin 9 r- . 
—7=— sm 6V + —v- ( v 2 cos 6>y - sin 6> y ) 
2V2/o 4/ 8 

cos6> . /?:•/, \ sin # 
— tv (cos 6» y + V2 sm V ) - cos Oy 

4/8 2V2/o 

cos 6* sin^ . /- „ ^ 

COS &y — — (cos b'y + V 2 sin 6V 

4 fo 



2v/2/o 



5fyp 7 (exp. 



(1.85 ±0.34) GeV" 1 
(1.31 ± 0.12) GeV' 1 
(0.60 ± 0.15) GeV 
(0.45 ± 0.06) GeV 
(0.70 ±0.03) GeV 
1.85 GeV" 1 
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-1 



-1 
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Table 2: Theoretical formulas and experimental values of the on-shell V-r], rj' electromag- 
netic vertex couplings defined in Eq. (|T2|). In Fig. ^| the theoretical values are plotted as 
functions of 9. We use 9y = 40.3° for the (f)-uj mixing angle. Experimental values are taken 

from IHH 



about 10 to 15%, as one may also infer from the cleanest channel uj — > ttj investigated in 
Sec. 2 (cf. Table p. A look at Fig. 0(b) shows that this is not the case for rj' — > pj. In that 
channel, however, we expect a contamination of the experimental value from non-resonant 
7T7T final states, which are produced by a completely different mechanism, which is related 
to the box-anomaly diagram (cf. e.g. ji~3f). For p — > 777 and uj — ► 777, Figs. ||(a) and (c), the 
experimental errors are too large to allow any serious restriction on 9. Within the error 
bars, however, our predictions agree with the experimental values, g^'-y, on the other 
hand, is very sensitive to the <p-uj mixing angle; for instance, reducing 9y by 1.5° brings 
the coupling within 20% agreement with the experimental data, i.e. just what we would 
expect. 

As seen from the figures, the most stringent constraint on 9 currently stems from the 
decay — > 777, which, as will turn out below, is only marginally consistent with the range of 
mixing angles found from other processes (remember however that we use 1 a experimental 
errors, and do not take into account theoretical uncertainties). The issue should be clarified 
when experimental errors will eventually shrink, and in particular we look forward to a 
measurement of the — > 7/7 mode, for which our prediction is only a factor 2 below the 
current experimental upper bound. 

5. Let us now turn to the investigation of other matrix elements of 77 and rj' that can be 
probed in our approach. We start with the light quark matrix elements, which determine 
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Figure 3: Theoretical (solid line) and experimental values (dashed-dotted lines, including 
experimental error) of the decay rate T(r] — > 3n ) as function of the pseudoscalar mixing 
angle. The short dashes give the theoretical error of the rate, which is dominated by the 
error of the quark mass ratio. 



the decay 77 — s- 37r. As in Ref. 

( I M75M 1 77 ) 
which yields the decay rate 

r(r^37r°) Vl > 



we findQ: 

1 /o cos 9 — v^/s sin 9 
~~ 75 v / 2/ocos# + f 8 sin 9 



( I S75S I 7] ) 



(21) 



4608vr 



3m 7r ) 2 5 r , 



m d - m u 



/|cos 2 g / /pcosfl- y/2/ 8 sinfl N 
cos + / 8 sine, 



(22) 

with a kinematical factor 5^ = 0.86. The most crucial ingredient in that formula is the 
ratio of quark masses r = (rrid — m u )/m s . There exist several possibilities to fix r from 
next-to-leading order calculations in chiral perturbation theory, as dicussed in 
One of them is based essentially on the analysis of 77 



14 



e.g. 

37r in Ref. \TB\. Since in that 
study 77 is necessarily treated as Goldstone boson, the value of r obtained that way is 
inconsistent with our analysis. For similar reasons, we also refrain from taking into account 
r determined from T(ip' — > J/if}T])/T(if}' — > J/ipir). We thus rely on the determination of r 
from pseudoscalar mass relations, which according to Ref. \TA\ yield r = 0.030 ± 0.005. 

In Fig. |3] we plot the decay rate as function of 9. Although the agreement with the 
experimental rate is very good for 9 ~ —17°, the theoretical errors are large and dominated 
by the uncertainty in the quark mass ratio. Actually the true theoretical error may be even 



larger and the rate be enhanced by final state interactions [U], [15 

In Ref. we had also considered the decay 77' — > 3tt. We have however since argued 
PI that this channel is considerably more complex and also receives contributions from the 
leading decay mode rj' — ► rjnn through rj-ir mixing. While this study is of great interest in 
itself and while a measurement of the charged mode 77' — > 7r + 7r - 7r° could shed considerable 



light on the mechanism underlying rj 
oiO. 



777T7T, it is of little help in determining the value 



6. We finally turn to the most important issue concerning the 77-77' system, namely the 
glue connection. Following the suggestion of Novikov et a!., Ref. [fT7 |, the glue matrix 



8 Eq. ( |2l| ) differs slightly from the corresponding Eq. (19) in Ref. B, where we had made the simplifi- 
cation / // 8 »l. 
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9 

Figure 4: Rj/^ = T(J/ip — > r]'j)/T(J/ip 777) as function of the mixing angle 9. Solid 
line: Rj/^p according to Eqs. ([17]) and (|23|) including experimental errors of /o an d /s (short 
dashes). Long-dashed line: the naive prediction (|2~4]). Dashed-dotted lines: experimental 
value with errors. 



elements are tested in the ratio of the decay rates T(J/ip — > r)(r]')j): 



m - r( j/v» - 777) 



( I 1 77' ) 



( 1 GG 1 77 ) 



™«/ m j/l, 



3 ' 



(23) 



Before turning to numerics, let us first comment shortly on the theoretical accuracy of the 
above equation. Actually the decay rates were calculated in exactly the same approxima- 
tion we used in deriving the gyp 1 given in the tables, namely dominance of the ground 
state and neglection of continuum contributions to the dispersion relations. It was also 
pointed out in Ref. [18| that the decay mechanism via the strong anomaly dominates only 
due to the smallness of the c quark mass and may be non-effective already in T decays, 
for which unfortunately no experimental data exist to date. It is for these reasons that we 
prefer to consider the ratio Rj/^ instead of the two decay rates separately: both radiative 
and continuum corrections are expected to cancel to some extent in Rj/^. 

In Fig. [| we plot Rj/^ according to Eqs. fllTD and ( p3|) and also the naive prediction 



R m = coM ) m " /m ^ 3 (24) 

as functions of the mixing angle. As experimental input we use R e f^ = 5.0 ± 0.6 Q. 
Our prediction ( |23|) is rather sensitive to the mixing angle. If we estimate the remaining 
theoretical uncertainty in (p3|) as ~ 10%, we obtain 9 = —(17.3 ± 1.3)°. It is clearly seen, 
on the other hand, that the naive curve is less consistent with the data and requires a very 
large mixing angle. To understand the departure of our formula from the naive prediction, 
we briefly evoke the elegant, but heavier formulation by Veneziano and collaborators ||. 
These authors deal with a simplified model, where 778-770 mixing is neglected. To apply their 
approach to the physical situation requires the introduction of three initial fields (instead 
of two in our approach): 

0o = 7^ (u-f 5 u + d~f 5 d + 575s) , 0g = -j= (uj 5 u + d-f 5 d - 2s7 5 s) , Q = ^ G Allv G Afxv . 

(25) 
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The definition of interpolating fields corresponding to physical particles requires the or- 
thonormalisation of these initial fields. While various procedures are possible (like in- 
troducing some "glue" Q in the definition of rf, as we did above when using the total 
divergences as interpolating fields), the choice of basis is essentially free, and the authors 
of Ref. |2j preferred to redefine the Q field, introducing G as belowQ This leads to (we 
follow closely the notation of 0): 

<f> v = b 8 cos 9<pl-b sin 9 <f%, </y = b 8 sin 9 <j>\ + b cos 9 G = Q - ■=^- </y - 

(26) 

where the coefficients hi and can be determined in terms of two-point functions. Ex- 
pressing the relevant Ward-Takahashi identities in terms of the effective action and taking 
two derivatives with respect to the gluon fields (g, with indices omitted below), we obtain 
for the (non- anomalous) octet axial current: 

T <t>lgg = b * sin ® T v'gg + b 8 cos 9T mg + XT Ggg = (27) 

with X a computable coefficient. Clearly, if the proper vertex between G and two gluons 
could be neglected (in the way the similar coupling between G and two photons is legiti- 
mately neglected in Ref. 0), we would recover the naive result (0). This step is however 
not permitted since the proper vertex between G (whatever its physical realization) and 
two gluons is precisely expected to be large. It should thus not be surprising that we have 
to depart from Eq. (p4|), which is difficult to reconcile with the data, as already noted by 



Gilman and Kauffman in Ref. []T0| . We will not dive deeper into this approach, except 
to mention that it is a clear generalization of ours (in our case the weight of glue in rj' is 
directly fixed by our definition of its interpolating field, instead as resulting from a further 
diagonalization between the three states). 

7. Having a rather comprehensive view of the rj and rj' system and its connection with glue, 
we now address other promising channels for future investigations, to wit D s decays. Here 
two processes compete to produce r](r]')X final states: The first one is the decay c — > s, 
leading to a ss pair, which hadronizes to rj or rj', while the second one can be described 
as cs annihilation into W and two gluons. It is suppressed by the Zweig-rule, but may in 
the present case gain importance through the large (0 | GG | 77(77') ) matrix elements. We 
should thus expect both an enhancement of the (77 + rf)X modes (unless the interference 
is destructive) and a large 77'/ 77 ratio, once the rates are corrected for phase-space and final 
state interactions. 

If we ignore for the moment final state interactions, the simplest way to evaluate a 
possible enhancement is by direct comparison to similar D decays. For instance the 
partial width for D — > K~p + (to which only the first diagram contributes) would, in 
absence of the glue contribution, be expected to be of the same order (after phase-space 
corrections) as the sum D s — > i]p + D s — > rj' p. Instead we find: 



T(D S - 7/p)(/v)- 3 + im - V p)(Pv)~ 3 (m Ds 



T(D -> K-p)(p K )-' 3 Wd , 



4.2 ±1.0, (21 



9 Note that the question if G is to be identified with any physical state, a glueball for instance, is left 
open in Ref. ||. 
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while for the branching ratios we find 



r(£>.->0») (pA _ 15±15 m 

which indicates both an enhancement of the decay rate and a large ratio for these important 
decay modes (of the order of 10% each). The situation is unfortunately less clear in the 7r 
associated modes, where the comparison with D now gives: 

T(D S -> 0r) (zy) -1 + T(D S -> ^(p,)- 1 = Lfi ± ^ 



T{D ^K-7T){p K ) 
T(D S - 



^ =3.0 ±1.1. (31) 



while 

t(d s -> 

More sophisticated analyses of non-leptonic Z) s decays were the subject of several stud- 
ies, e.g. |fn% EH, PHI, which, however, all relied on the naive picture of 77-77' mixing without 




taking into account the anomaly. In Ref. |20| it was pointed out that in order to account for 



the data the non-leptonic rates involving 7/ need to be enhanced by an extra factor of up 
to three when compared with the corresponding D decay amplitudes, even after including 
final state interactions. And the global fit done in Ref. |21| for all non-leptonic D and D s 



decays including final state interactions, W exchange and annihilation contributions fails to 
explain the large branching ratios of D s into tj'tt and rj'p. The enhancement needed is how- 
ever well in line with our expectations from the extra Zweig-suppressed (gluon-mediated) 
amplitudes. 

Very recently, the CLEO collaboration has measured the theoretically cleanest semilep- 
tonic channel and finds |22j 

T(D -> rj'eu) = 

T(D s ^ V eu) 1 ; 

Defining the relevant formfactors as 

( V I s7m c \ D s) = f+(q 2 ) (pd s + P v )^ + f-(q 2 ) (pd s - p v )p (33) 
with q = pd s — p v and analogously f± , CLEO's measurement yields 

±±y- = 1.14 ±0.17 ±0.13, (34) 

where we assumed monopole form factors^ with a pole at q 2 = m 2 D * . The naive mixing 
model, on the other hand, predicts 

fl(q 2 ) _ cos 9- vising 

fl(q 2 ) sin^ + v^cos^' 1 J 

which translates into 9 = — (13.5±4.7)° and is only marginally consistent with the deter- 
mination from 77 — > 27 



10 We do not expect the anomalous contribution to change drastically the q 2 dependence of the form 
factors, which in D — > Kev is a monopole experimentally. 
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As a temporary conclusion to this section, we would like to add that, while this sector 
is clearly difficult to investigate both theoretically and experimentally, it could prove very 
fruitful for our understanding of the role of anomalies and the importance of gluons. There 
is also the strong suggestion that D s decays, not unlike the J / ip radiative decays, constitute 
a "glue- rich" channel, where search for possible glueball states should be considered. We 
hope to come back to that issue in future work. 

Let us finally remark that apart from understanding the structure of the rj and rf 
particles, the size of their connection to glue is of importance to investigate possible glueball 
candidates [p3fl . In particular, the candidates /o(1500) from crystal barrel |24| and /o(1590) 



ZS| are seen to decay more frequently into rji]' than into r/r] (after phase-space correction), 
a ratio which we expect to be given by 

r(/o -» v'v) 



r(/o -> vv) 



(0\GG\rf) 



( 1 GG | rj ) 



CM 

(36) 



Pryrf 



n CM 



where P^fn are the respective centre of mass momenta. This relation is however highly 
momentum dependent, since we are close to the threshold, and obviously needs to be 
integrated over the particle width. 

In conclusion, the above study gives a reasonably consistent description of the rj-rf 
system for values of the mixing angle 9 between —20 and — 17 degrees, taking into account 
anomalies, as tested from various channels, respectively depicting eletromagnetic proper- 
ties, light quark and gluon content. Our approach will be tested by a general improvement 
of the experimental measurements, including the V — > P7 and P — > V7 rates, and in par- 
ticular by the observation of <fi — > 7/7. The investigation of "glue-rich" channels, namely 
the traditional radiative decays of J/ip, but also, as advocated here, the D s decays, will 
also yield a better understanding of this system and at the same time help in the search 
for glueballs, which might solve the puzzle. 
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